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Abstract 

A comprehensive review is presented of the progress made in further developing the ghost-free Analytic Approach 
to low-energy QCD since "QCD-97" meeting. It is now formulated as a logically closed "Analytic Perturbation 
Theory" algorithm. Its most essential feature is nonpower functional expansions for QCD observables. 



1. Intro 

The claim is that current practice to use 
expansions for QCD observables in powers of 
the same effective coupling function a s , i.e., 
[a s (Q 2 )] k , [a s (s)] fc or [a s {r)] k , in different pic- 
tures (Q 2 , s , r ...) as, e.g., it is implicitly recom- 
mended by PDG[1], is neither based theoretically, 
nor adequate practically to low-energy QCD. 

Instead, one should use non-power function- 
al expansions with non-power sets of functions 
{Ak(Q 2 )} ,{2lfe(s)} and {H fc } , related by some 
suitable integral transformations. 

To comment this, consider a class of homoge- 
neous Linear Integral Transformations (LITs) 



fE -> fi(Xi) 



dQ 2 
Q 2 



(1) 



that describe transitions from the Euclidean pic- 
ture to Minkowskian and Distance ones: 

Je(Q 2 ) - f M (s) ■ f E (Q 2 ) - MiA 2 ) . 



Email address: shirkovdOthsunl.jinr.ru (Dmitry 
Shirkov). 



Being applied to renormalisation-invariant func- 
tions they could relate, in particular, effective cou- 
plings properly defined in these pictures 



ue{Q 2 ) -> a M (s) ; asiQ 2 ) —> a> D (r) . (2) 



Meanwhile, generally, LITs (I) (and reverse 
transformations) can be used only within a class 
of functions f t that are free of nonintegrable 
singularities. For example, a common expression 
for QCD effective coupling a s with powers of 
L = ln(Q 2 /A 2 ) in denominators (like eq.(9.5) in 
ref.[l]) is not suitable here. Due to this, in our 
further discussion we use some adequate models 
for a s , free of unphysical singularities. 

Moreover, LITs change the form of observable 
expansion; being applied to powers of some a{Q 2 ) , 
they yield expressions that are not powers of the 
a{Q 2 ) image. They transform expansions in cou- 
pling powers a k into expansions over non-power 
sets |A(Q 2 )} ,{2tfc(s)} and {M^ 2 } ■ 

We give explicit examples of these sets and con- 
struct non-power expansions for observables that 
demonstrate quicker convergence and reduced 
rcnorm-scheme sensitivity. 
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2. Different pictures in HEP 

2.1. Various HEP representations 

An essential property of any successive theoret- 
ical scheme for physical phenomena is the possi- 
bility to use various pictures for their description. 
In HEP, these are space-time, energy-momentum, 
eikonal and some others. Below, we discuss the 
momentum-transfer (Euclidean part of the energy- 
momentum 4-space defined by Q 2 = Q 2 — Q 2 , — 0)' 
the ccntcr-of-mass-energy squared (Minkowskian 
part of the energy-momentum 4-space with s = 
<7q — q 2 > ,) as well as the distance (space-like 
part of space-time r 2 = x 2 — x\ > ) pictures. 

Remind, first, that Euclidean momentum- 
transfer picture was, historically, a natural 
scene for formulating [2a-c] (see, also review 
Chapter "Renormalization Group" in [3]) basic 
renormalisation-group (RG) relations and equa- 
tions. Here, an important feature is that scalar 
propagator amplitudes depending on a real argu- 
ment Q 2 > are real functions, and their particu- 
lar values can be expressed via real parameters Zi 
of finite Dyson renormalisation transformations 1 . 
Just due to this, we use the Euclidean picture with 
Q 2 > as a launch-platform for LITs (1). 

2.2. Pictures related by integral transformations 

Transition from an Euclidean function Je{Q 2 ) 
to the Distance one foil* 2 ) can be performed with 
an appropriate Fourier transformation. For exam- 
ple, in the three-dimensional case one has 

fE(Q 2 )^¥[f E }(r 2 )^f D (r 2 ); 



¥[f E ](r 2 ) =r dQ sin(Qr) f E (Q 2 ) . (3) 
Jo 

At the same time, in current QCD analysis, 
people use a relation between the Euclidean and 



1 Analogous possibility is provided by the Euclidean space- 
time picture first used by Dirac [4] for introducing space 
distribution e(r) of the electron charge. 



Minkowskian pictures based on definition of the 
Adler function 

We treat this relation as integral transformation 

R ( s ) -> D (Q 2 ) = D [R] (Q 2 ) , (4) 

that transforms a real function R(s) of the 
Minkowskian domain (real s > ) to a real func- 
tion D(Q 2 ) of the Euclidean one (real Q 2 > 0). 
The reverse transformation M = has the form 
of a contour integral. 

Both the transformations (3) and (4) are consis- 
tent with RG-invariancy and are non-consistent 
with existence of unphysical pole at A 2 . 

3. Ghost-free APT construction 

3.1. Ghost-free QCD coupling in Minkowskian 

In the early 80s, by summing the 7r 2 -terms aris- 
ing in the course of logs branching 

'Q 2 ' 



1= In 



A 2 



^ (£) - - 



ITT 



under the transition Q 2 — > — s Radiushkin, Kras- 
nikov and A. Pivovarov obtained remarkable ex- 
pressions for Minkowskian QCD coupling and its 
"effective powers" . This operation , in the 

one-loop case for — l//3o£, results in the ex- 
plicit form [5]: 



We \ 

a M \ s ) — ~q arccos 



1 



~2 ' 



(5) 



1 7T 

= — — arctan — 
l>o Pq-k L 



and for square and cube of ct-p in expressions [6] 
1 ..m, . 1 



(3 2 £ 2 



2l W ( S ) 



(i 2 [L 2 + ^] 
L 



(6) 



P 3 o \L 2 



which are not powers of ofa ( s ) ■ As it was noticed 
later [7] 2 , they are ghost-free. 



See also Rcf . [8] . 
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3.2. Ghost-free QCD coupling at Euclidean 



In the Euclidean momentum-transfer domain 
analogous ghost-free construction for QCD cou- 
pling and observables was proposed [9] in the 
mid-90s. This Analytic approach was based on 
imperative of analyticity in the cut Q 2 -plane. Its 
algorithm, first realized [10] in QED, combines 
RG invariancc the Kallen Lehmann spectral rep- 
resentation with — see, Sections 53 and 54 in [3] 
and the first of eqs.(9) below. 

In the simplest case of 1-loop QCD coupling, 
the analyticization operation A yields analytic Eu- 
clidean QCD coupling 



>/3o4 1} (Q 2 ) 



A 2 



Q 2 - A 2 



In 



91 
A 2 



Analogous expressions with "subtracted" high- 
order pole singularities are valid for powers of afp ; 
e.g., for coupling squared 



1 

T 



i 



Q 2 A 2 



(Q 2 -A 2 r 



(7) 



By construction, new analytic functions 
{A k } = Ai{= a E ),A 2 ,- : 
■jt are free of unphysical singularities; 
if are finite in the IR limit with 

a B (0) = l/A), A> 2 (0) = 0; 

it include non-perturbative structures; 

it in the weak-coupling UV limit — > a* ~ [(]~ k ■ 

These properties remain valid for higher-loop 
cases. In particular, IR limiting values (7) pro- 
vide the basis for remarkable stability of the Ak 
functions behavior in the low-energy domain. On 
Fig.l, one can see that while the 2-loop curve 3 
for oie(Q 2 ) differs from the 1-loop one in the few 
GeV region by 10-20%%, the 3-loop approxima- 
tion practically coincides with the 2-loop curve, the 
difference being of the order of 2-3 per cent. The 
same is true for Minkowskian coupling. This, in 
turn, yields low sensitivity with respect to change 
of renorm-scheme (RS). 



3 For analytic expressions and numerical tables of 2-loop 
functions see papers [9,lla-e]. 
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Fig. 1. Space-like and time-like global analytic cou- 
plings in a few GeV domain with A^= 3 ) = 350 MeV . 
Here, ctE,M In 1-loop case are drown by 
dash- and- dotted lines and for 2- and 3- loop cases by 
solid curves. 

3.3. Minkowskian vs Euclidean; the APT 



As it was noticed later [12], both the sets of Eu- 
clidean and Minkowskian functions are connected 
by LITs A k (Q 2 ) = D[a fc ] , 2t fe (s) = R [Ak] ■ In 
the 1-loop case [7] 

1 L 1 1 A 2 

— arccos = 

7T Vi 2 + 7T 2 



1 



+ ■ 



ln(Q 2 /A 2 ) 
Q 2 A 2 



Q 2 - A 2 



1 



L 2 + 7T 2 ' 



_ln 2 (Q 2 /A 2 ) (Q 2 -A 2 ) 2 _ 

Accordingly, for an Euclidean RG-invariant func- 
tion Ie(Q 2 ) ~ initially expandable in powers of 
the coupling constant a s - there appears a rela- 
tion with its Minkowskian image Im{s) — R[-/s] 



'M 



i E (Q 2 )^J2 l "MQ 2 )- 

k 

(8) 

Thus, starting with common RG-improved per- 
turbation expansion for an Euclidean function 
Ipt(Q 2 ,ct s ) , after applying operation of analyti- 
cization A , and of tt 2 -summation J2ir 2 we ar ~ 
rive at non-power functional expansions related 
by (8) with numerical coefficients %k ■ The whole 
construction has been named [7,12,13] the An- 
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alytic Perturbation Theory (APT). It was duly 
elaborated [12] for real QCD with various numbers 
/ of active quarks. The logic of its structure is 
presented in Fig. 2. 



RG + PT 




«M(s) 




a E {Q 2 ) 


F 


a D (r) 






A k (Q 2 ) 




H(r" 2 ) 


K 




Minkow 


Euclid 




Distance 



Fig.2 Logic of APT scheme. 



3.4. Sketch of the global APT algorithm 

The most elegant form of the APT formalism 
uses a spectral density p{a) taken from pertur- 
bative input. Then all involved functions in the 
above-mentioned pictures look like 



oo 

1 f do- 

- / — Pk(o-), 
ir J a 



(9) 





In the 1-loop case, 

^ 1) = M^J ; L * = ln G&) (10) 

with higher spectral functions expressed via 
by simple iterative relation 



dL„ 



(11) 



that implies analogous iterative relations for Ak 
and 2lfe . As it has been noted above, eqs.(10) and 
(11) were generalized for a higher- loops case with 
transitions across heavy quark thresholds and were 
successively used [13] for fitting various data. 



Let us emphasize an important ansatz built- 
in the presented APT construction: spectral 
functions, like (10), are defined on the basis of 
RG-improved perturbation theory. Due to this, 
they contain no additional parameters and non- 
perturbative elements. Hence, this construction 
can be considered as a "minimal APT scheme" . 

Here, it can be added that some detail of the 
APT construction in the distance picture was re- 
cently considered in our paper Ref.[14]. 

4. Non-power expansions for observables 

To summarize APT for observables, we state 
that instead of usual expansions in powers of the 
same universal QCD effective a s function 

{a k s (Q 2 ,a s )}, {a k s (s,a s )}, {a k s (l/r 2 ,a s )} 

one should use "perturbatively motivated" expan- 
sions in terms of non-power sets 

{A k (Q 2 ,a s )} , {K k (s,a s ) = R{A k }(s)} , 

|H fc (l/r 2 ,a s )=F[A](l/r 2 )} 
with some of them being non-analytic at a s = . 

4.1. Properties of non-power functions 

Review now the main properties of nonpower 
sets {*4fc},{2lfc} . Qualitatively, all three 

are similar: 

I. They consist of ghost-free functions 

II. First functions, new couplings, oe, ojm, 
<^> are monotonic; 

^> In the IR limit, they are finite = l//3n — 1.4 ; 
In the UV limit <~ 1/ In x ~ a s (x) . 

III. All other functions (k > 2) : 

9 start from zero A k (0), » fc (0), N fe (0) = 0; 

V in the UV limit - l/(lna;) fc - a k s {x) . 

•v 1 2nd ones, Ai , 2I2 , ^2 obey max at <~ A 2 . 

•s? Higher ones k > 3 oscillate at ~ A 2 and obey 

k — 1 zeroes 4 . 

Just two last properties play an essential role is 
the RS-inscnsitivity and in better convergence of 
expansions in the low-energy region. 



On mathematical nature of these sets see papers [15,16]. 
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Fig. 3. "Distorted mirror symmetry" for global expan- 
sion functions {A k (Q 2 )} and {2l fc (s)}. All the curves 
here correspond to exact two-loop solutions expressed in 
terms of Lambert function. One can see that 3rd func- 
tions arc oscillating in the low-energy region. There, they 
are much less than third powers (dot-and-dash lines) of 
related coupling functions. 

It is curious that the mechanism of liberation of 
singularities is quite different: for Ak , in the space- 
like domain, it involves subtraction of nonpertur- 
bative, power in Q 2 , structures. Meanwhile, in the 
time-like region, it is based only on resummation 
of "7r 2 - terms" . Graphically, nonperturbative ana- 
lyticization in the Q 2 -channel looks like a slightly 
distorted reflection (as Q 2 — ► s = — Q 2 ) of the per- 
fectly perturbative "ir 2 -summation" result in the 
s-channel. The effect of "distorting mirror" (first 
discussed by Milton-Solovtsov [7]) is illustrated in 
Figs. 1 and 3. 

4.2. Better convergence and RS insensitivity 

Due to a quite different from the usual oscillation 
behavior of higher (with k > 3 ) functions, the 
APT nonpower expansions 

cU(Q 2 ) = di a E {Q 2 ) + d 2 A 2 (Q 2 ) + d 3 A 3 + ... 
r w (s) = dia M {s) + d 2 s &2(s) + d 3 2t 3 (s) + ... 

obey better convergence in the hadronization re- 
gion. We demonstrate this by a Table from paper 
[13]. All results are given in the MS scheme. 



As it follows from the comparison between the 
3rd columns for usual (PT) case and for the APT 
one, the relative contribution of the 3rd ( ~ ) 
term drastically diminishes; it is not any more of 
the same order of magnitude as the 2nd term. 

Table: Contributions in %% of 1- , 2- 



and 3-loop terms into observables 



Process 


PT 


APT 


1st 


2nd 


3rd 


1st 


2nd 


3rd 


GLS Sum Rules 


65 


24 


11 


75 


21 


4 


Bjorkcn Sum Rules 


55 


26 


19 


80 


19 


1 


Inch V. r-decay 


55 


29 


16 


88 


11 


1 


lOGcV e+e- -> hadr 


96 


8 


-4 


92 


7 


.5 


Z — ► hadrons 


98.6 


3.7 


-2.3 


96.9 


3.5 


-.4 



In turn, this leads to essential reducing of 
dependence of the results on change of the 
renormalisation-scheme prescription. For further 
numerical examples of reduced RS sensitivity see 
papers [17,18,19]. One more fresh illustration is 
provided by recent analysis of the pion form factor 
[20] presented by N.Stefanis at this meeting. 

4.3. Issue of the APT non-uniqueness 

The described minimal APT scheme obeys one 
specific feature: it does not contain any new pa- 
rameters. Like in usual QCD case, there is only one 
adjustable parameter; it can be chosen as h(f=5) 
or a s (mz ) • Aesthetically, this feature looks very 
nice. In particular, as a thorough numerical anal- 
ysis reveals, within the APT framework with its 
two different coupling functions % and a_g , it 
turns out to be possible to correlate a great bulk of 
QCD data more successively [13] (with smaller x 2 
value) than by a standard algorithm with the only 
a s - see the paper [21]. 

Meanwhile, there exist nonperturbative evi- 
dence that in the 1 GeV region QCD coupling 
could reach values close to 0.6-0.9, [22a-c] (see also 
our review [23]) which is impossible for the mini- 
mal APT effective couplings at realistic A values. 

In APT, there are several possibilities for insert- 
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ing adjustable "low-energy parameters" . The sim- 
plest - from physical point of view - one, is to at- 
tribute some effective low-energy mass to partons. 
This idea is not a new one — see, e.g., refs.[24a-c]. 
In the APT context this trick was successively used 
by Baldicchi and Prosperi [25] in calculating mass 
spectrum of light mesons. 

Another way is to insert some non-perturbative 
structure into spectral density p(a) in eqs.(9). An 
example of this kind can be found in Rcfs.[26]. 
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